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Abstract. The self-organized hydrodynamic models can be derived from the kinetic version 
of the Vicsek model. The formal derivations and local well-posedness of the macroscopic 
equations are done by Degond and his collaborators. In this paper, we rigorously justify this 
hydrodynamic limit. 
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1. Introduction 

It has been an active research area in recent years on the modeling of interacting agent 
systems arising in nature, such as bird flocks, fish schools, herds of mammals, etc. They 
provide fascinating examples of self-organized systems which are able to produce large scale 
stable coherent structures. Among these models, the Vicsek model [12] has received particular 
attention due to its simplicity and the universality of its qualitative features. This model 
is a discrete particle model which consists of a time-discretized set of ordinary differential 
equations for the particle position and velocities. The velocities are assumed to be of constant 
norm and updated according to an alignment rule, i.e. each agent tries to align its velocity 
to that of its neighbors in some sensing region. Some angular noise is added to account for 
stochastic fluctuations. A time-continuous version of this model and its kinetic formulation 
are available in Oi. From the the time-continuous Vicsek model to this kinetic model is 
rigorously justified in [1]. 

In [8], Degond-Motsch derived formally the hydrodynamic limit of the time-continuous 
Vicsek model through an asymptotic analysis of its kinetic version [T]. The resulting model is 
a system of balanced equations for the density and mean velocity orientation (or polarization 
vector). This model was later called the Self-Organized Hydrodynamics (SOH). Furthermore, 
in [7], Degond-Liu-Motsch-Panferov derived the modifications of the SOH model by the in¬ 
troduction of the attraction-repulsion force, different scaling assumptions about the size of 
the sensing region which involve a higher level of non-locality. As proved in |8], the strict 
combination of alignment and noise results in the appearance of a pressure term in the SOH 
model. Introducing an attraction-repulsion force and spanning various scaling assumptions 
on this force and on the size of the sensing region, they obtain in [7] a variety of effects that 
are not encompassed in [8]. 

In [7], besides the derivation of the macroscopic equations, they prove a local-in-time exis¬ 
tence theorem in the 2D case for the viscous model (when the non-local effects are retained) 
and in the 3D case for the inviscid case (when the non-local effects are omitted). Both proofs 
are based on a suitable symmetrization of the system and on the energy method. 

The main concern of the current paper is on the rigorous justification of the hydrodynamic 
limit from the Self-Organized Kinetic (SOK) system to the Self-Organized Hydrodynamics 
(SOH). The main challenge is the lack of conservation laws. To overcome this difficulty, the 
Generalized Collision Invariants (GCI) are employed to derive the macroscopic equations. We 
start from the solution constructed in [7] and proved that there exists a class of the solutions 
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of SOK uniformly on a time interval on which the solutions of the SOH are constructed, such 
that their hydrodynamic limits are the SOH. We employ the methods developed in the fluid 
limits of the Boltzmann equations, such as [2] and m- 

The paper is organized as follows. In section 2, we introduce the self-organized models: 
the time-continuous Vicsek model, the kinetic model and the formal hydrodynamic limits. 
The main results are stated in section 3. Some mathematical preliminaries are introduced 
in section 4. In section 5, we estimate the term appear in the expansion, and the main 
estimates for the remainder equation is presented in section 6. In the last section, building 
on the estimates in the previous two sections, the proof of the main theorem is completed. 

2. The Self-organized model 

2.1. Self-propelled particles interacting through alignment. Let Xk{t) G and 
Vk{t) G be the position and velocity of the A:-th particle at time t. The time-continuous 
version of the Vicsek model is written as follow. 

(2.1) Xk{t) = cXk{t), 

(2.2) dVk{t) = Py± o {vVk{t) dt + VwdB^) , 

(2.3) = ^ y., 

' h\Xj-Xk\<R 

where for V G Pyx = Id — V(^V is the orthogonal projection onto the plane orthogonal 

to V. The equation (|2.2I1 takes the form of a stochastic differential equation (SDE). The 
projection operator Pyx ensures that the resulting solution Vk(t) stays on the unit sphere, 
provided that the SDE is taken in the Stratonovich sense (which is indicated by the symbol 
o). The hrst term inside the bracket is the interaction. It corresponds to a force acting in 
the direction I 4 of intensity The second term is a white noise consisting of independent 
Brownian motion in of intensity y/2D. Then, let V ^ 00, a mean field model is 
obtained. This model is described as follows. 

2.2. Mean-field model and scaling. The mean-held model describes the evolution of the 
one-particle distribution function /(t, x, v) at position x G R"", with velocity v G at time 
t >0. The model is written as: 

(2.4) dtf + cv • VJ = -V, • {Fff) + DAJ , 

(2.5) Ff{t,x) = uP^±Vf{t,x), i;/(t,x) = A!A^ 

\Jf[t,X)\ 

(2.6) Jf{t,x)=[ P:((^)/(t,y,a;)wdyda;, 

where the constants c, v, D are the same as in (12.ip and (j2.2p . The equation (j2.4p is a Eokker- 
Planck equation. The left-hand side expresses the rate of change of / due to the spatial 
transport of the particle with velocity cv while the hrst term at the right-hand side denotes 
the transport in velocity space due to the interaction force Ff. The last term at the right- 
hand side is a velocity diffusion term which arises as a consequence of the Brownian noise in 
particle velocities. Note that since v lies on the sphere. Ay stands for the Laplace-Beltrami 
operator on the sphere. The derivation of the mean-held model ()2.4I) . ()2.5p and ()2.6p from 
the discrete system (HTP, (I22|) and ([231) has been justihed in [I]. 

After nondimensionalization process, (for details, see [6]), we can write down the scaled 
self-organized kinetic (SOK) model, 

dtf + V- V,/ + % Vy • (PyXP//) = i (-Vy • (PyXDj/) + dAy/) , 


(2.7) 
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where the parameter e denotes the mean free path, i.e., the distance needed by a particle to 
make a finite change in direction of motion due to the interaction force, and r/o = 0,1 indicate 
the inviscid and viscous case, respectively, and d > 0 is a constant. Furthermore, the local 
current density jf, the local mean direction flj-, and the quantity Dj are defined respectively 


as follows: 



(2.8) 


= [ vfdv, 


JS"-1 

(2.9) 

^f{t, x) 

del jf p;n-l 

~W\ ^ ’ 

(2.10) 

Df{t,x) 

X p ^^jf 

^ li/I ’ 


2.3. Basic Properties of the SOK Model. In this subsection, we list some basic prop¬ 
erties of the self-organized kinetic (SOK) model, and we refer the readers to [U O |9] and 
references therein for more details. 

First we propose a hypothesis: = 0, if = 0. Define the collision operator Q by 

(2.11) Q(/) = -V^-(P,xDj/) + dA,/. 

We mention that Q{f) = dA^f in the case \ jf\ = 0 by the above hypothesis. 

Next we describe the equilibria of Q, which are expressed by the von Mises-Fisher (VMF) 
distributions with respect to the local mean orientation D € namely, 

f ={/| Qif) = 0} = {v^ pMniv), Vp € M+, D G 


where the VMF distribution is defined as 

(2.12) Mn{v) = Z-^ exp(^^^) 

with a constant Zd = exp(^)du independent of fi. The VMF distribution enjoys 

the following properties: 

i) Mq(v) is a probability density (i.e., /^ggn-i Mq(v) du = 1); 

ii) The first moment of Mq{v) satisfies 



1 


QMq{v) du 


ciD, 


ci(d) 


IgS^-i(^-^)exp(^)du 

/„esn-iexp(^)du 


where the coefficient ci G [0,1] denotes the order parameter in the study of phase 
transitions. 

Note that the formula = V.y(u • D) ensures that the collision operator Q can be 

rewritten as 


Q(/) = d^. 


which results in a dissipation relation 




f 


Mn, 


(2.13) 




2 

du < 0. 


This implies that Q{f) = 0 is equivalent to / G T. 

One of the main difficulties to derive the macroscopic equations of the SOK model is that 
it obeys only the conservation law of mass. To recover the missing momentum conservation 
related to the quantity Degond-Motsch introduce the concept of the “Generalized 

Collision Invariants” (GCI) in [8], see also [Umil] its applications. 










4 


N. JIANG, L.-J. XIONG, AND T.-F. ZHANG 


Definition 2.1 ([8]). For any given 12 E §” the linearized collision operator Cq, is defined 
as 


(2.14) 


Cnf = A,/ 


iv, • (P,xl!/) = V, 




The Generalized Collision Invariants (GCI) are the elements in the null space oj Tq,: 
(2.15) M{Cn) f Cnf fidv = 0, V/ such that 12^ = ill 

= {fi\ Cnfi{v) = A - V with ^4 • 12 = 0} 

= {u I—>■ h{v ■ 12)^4 ■ V + C with C E M, A E M”, and A • 12 = 0} , 


where the operator is the adjoint of the linearized operator Cn, which takes the form, 
C*nfi = - ^12 • Vyfi = ■ {MnVyfi). 

Here, h{v ■ 12) = h{cos9) = with g being the unique solution of the elliptic equation 
Cf^g{9) = sin0 in the space V, where 


^h9{S) = -sm‘^ "6*6 "“"^(sin'^ 9e''°d%' ( 0 )) + !!-^g{9), 

at) sin 9 

y = 1^1 (n — 2) sin2 E L^(0, vr), sin2 6*5 E Pg (0,7r)| . 


Using the GCI, the macroscopic equations of 12 can be derived, as stated in the following 
subsection. 


2.4. Formal derivation of the self-organized hydrodynamics. To study the macro¬ 
scopic limit of the self-organized kinetic (SOK) model, we rewrite (12. 7p as follows, 

(2.16) dtf + V ■ V,f + r?oV, • {P.rDfef) = ^Q{f). 

We seek for a special class of the solutions of (j2.16p of the form: 

(2.17) /*^ =/o + ^/i + e^/l) 
with the restriction: 

( 2 . 18 ) jfe = jfo . 

As a consequence, 12je = 12/g = 12o • Using this llg, we can define the Generalized Collisional 
Invariants as in Lemma EH Under the restriction (I2.18p . the nonlinear equation (12.161) 
becomes linear as follows. 

(2.19) dtr + V • + % V, • {P^vDof) = -^Cn,r , 

where Dg Df^. We plug (I2.17P into the equation (I2.19D . and collect the same orders, which 

gives: 

Order 0(j) : The leading order is 

( 2 . 20 ) o = dCnofo = Q{fo)- 

Recalling the equivalence between Q{f) = 0 and / E T, the equation (|2.20p implies 

(2.21) fo{t,x,v) = po{t,x)Mnf,{t,x){v), 

for some function pQ = po{t, x). In the rest of the paper, we use the notation Mg = ■ 

Order 0{1) To determine the equations satisfied by the macroscopic variables (pg,12g), we 
consider the order 0(1): 
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We require the part of fi in vanishes, and can solve the part in , 

(2.23) h = {dtfo + V • V,/o + r?o V, • {P,±Dofo )), 

under the solvability condition that the right-hand side of (I2.22p lies in (for details, 

see section 5). This gives the following macroscopic equations satisfied by (poj^^o): 


(SOH) 


' dtp + ciVx ■ (pf2) = 0, 

< p{dt^ + C 2 O • Vx^) + dPQ±VxP = C3Pq±Ax{p^), 

.M = i, 


with the coefficients 

_ cos 0 exp(^^) sin^-2 0 d0 

“ /(;exp(^)sin"-2 0d0 ~ 

C 3 = Pok[{n - l)d + C 2 ]. 


£ cos eh{cos e) exp(^^) sin” 9 d0 
£ /i(cos 9) exp(^) sin” 9 dO 


This is the hydrodynamic model, which we call self-organized hydrodynamic (SOH) system. 
We refer to [a Eli for the derivation of the SOH system and omit the details here. Now the 
equation of the remainder /| is 


dtfi + ^ + Vo • {Py-'-k^ofi) ~ 

= - {dtfi + V • Vxfi + r?o V, • {P,±Dofi)} . 


For notational simplicity, we set d = 1 in the rest of the paper. By setting /I 
equation (I2.24|] is reduced to the equation of /|: 


(2.25) 
where 

(2.26) 


dtfl + V ■ Vxfl + % V, • P,±Dof, 


+ -To/I = hofi + -hi 

e e 


A/ = -^Too(^o/) = • (MoVyf) 


/IMo, the 


(2.27) 
and 

(2.28) 


^ + ho V, • (P„xZ)o/i)] , 

Mq 

ho = —— [dtMo -I -V ■ VxMo + vo'^v- {Py^DoMo)] 
Mo 

= -[v ■ dt^o + V ®v : + VoDo ■ Py^^o] ■ 


In the rest of the paper, we work on the remainder equation (|2.25l) . 


3. Main results 


In this section, we state our main result. We first introduce the existence result of Degond- 
Liu-Motsch-Panferov, on which our result is built. First, we introduce the Cauchy problem 
of self-organized hydrodynamic (SOH) system: 


(SOH) 


dtp + ciVx ■ (pll) = 0, X E T, t > 0, 

p{dtVt -h C2O • Vj,!!) + dP^±VxP = csP^^Axipil.), 

\n\ = i, 

, p\t=o = p'^ > CO > 0, n\t=o = fi'”, |fi'”| = 1 . 


Here T denotes or T^. 

On the one hand, the SOH system evidently bears many similarities with the isentropic 
compressible Navier-Stokes (NS) system. And on the other hand, it also has some differ¬ 
ent properties. The hrst important difference is that the SOH system obeys the geometric 
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constraint |r2| = 1 which requires the velocity O to be of unit norm. The second important 
difference is that, generally speaking, the coefficients c* (z = 1,2) are different. Indeed, we 
have 0 < C 2 < Cl (the equality holds iff d = 0 ), see [6l [8] . 

The existence and uniqueness of the solution to the Cauchy problem of SOH system have 
been established in [7], we quote the results as follows. 

Theorem 3.1 ([7]). Let T = T” , given > cq > 0, n*” E i.e. 

• for n = 2, n®"’ = (cos , sin E [0, 27r]; 

• for n = 3, n®” = (sin 0®” cos sin d®*® sin , cos d®”), 0®®® E [0, vr], f/?®®® E [0, 27r]. 
Furthermore, assume sind®®® > 0 , and C 3 = 0 , (i.e. r]Q = 0 .) 

with (p®*®, (^®’®, 0®®®) E Id®®®, m > f+ 1 ■ Then, there exists T > 0, such that the Cauchy problem 
of SOH system with initial data (p®®®,!!®®®) has a unique solution (pojdlo) € T°°([0, T], dd®®®(T))n 
ddi([ 0 ,T],dd®'®-HT)), andpo> 0 . 


Our hydrodynamic limit result builds on the above theorem, so make the same assumptions 
on the initial data, namely, 

Assumption (A): 

(1) In 2D case, we consider C 3 > 0 including both viscous and inviscid cases. The initial 
data (/ 9 ®®®,D®®®) E dd®®®(T) are smooth enough as required. Besides, p®®® has a positive 
low bound. 

(2) In 3D case, we consider only the inviscid case C 3 = 0, which arises from the coefficient 
rjQ = 0. (Indeed, C 3 = rjtfc^). Besides, p®®® has a positive low bound and D®®® 7 ^ (0,0,1) 
(corresponding to the above constraint sind®®® > 0 ). 

Now we state the main result of this paper on the hydrodynamic limit from the self- 
organized kinetic (SOK) equation to the self-organized hydrodynamic (SOH). 


Theorem 3.2. Let (po,^^o) G L°°([0,T],H'®®®(T)) n ^^([O,T],H'®®®-i(T)) (m > 13) he the 
solutions provided by Theorem \3.1\ to the Cauchy problem of the self-organized hydrodynamic 
system (jSOHp with initial datum (p®®®,D®®®) satisfying the assumption (A). Let /o = pqMq^, 
and fi{t,x,v) E Af'^{Cng) be the unique solution given in (j2.23p . 

Furthermore, assume /^’®®®(x,z;) = p®®®(x)MQin( 2 ,)(u) + sfi{0, x,v) + u) with the 

bound ||/|’*’"(x,n) 11^12 ^ < C. 

Then there exists eo > 0 such that, for any e E (0,eo)? l-he self-organized kinetic equation 
(I2.16P admits a unique solution f‘^{t,x,v) & C{[0,T]-, H'^{T x§^~^)) in the class jfe = j, of 
the form 

(3.1) f{t,x,v) = po{t,x)Mn^(^t^,^){v) + efi{t,x,v) -\- e‘^fl{t,x,v), 


with f 2 {t,x,v) satisfying 

(3-2) +£||/|IIhi,„ 


3 

62 


<C, 


where the constant C is independent of e ^ (0,eo) o,nd t E [0,T]. 


Remark 3.1. Indeed, we can get a more generic result about the higher order diffusion 
expansion, 

(3-3) /^ = PO-^Qo + + ■ ■ ■ + ^” ^fn-l+A^fni 

if it holds initially with the bound \\fn^^\\H^„ < C. For i = 1,2, - ■ ■ ,n — 1, fi are determined 
by the equation 


Tnofi = dtfi-i -\- v ■ Vxfi-i + VoVy ■ {Py±Dofi-i), 


and the n-order remainder ff satisfies 

(31) E 




Iff5 


<C, 


0<s<n 
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with the constant C independent of e G (0,eo) OTirf t E [0, T]. 

Proving this high order result works exactly the same as for Theorem \3.^ and relies heavily 
on a high order a priori estimate (I6.2ip . see the discussions in Remarks \6A\ and \5.1\ below. 


4. Preliminaries 

In the following context, we will do estimates in the weighted spaces L'^ {Mq dv dx). For 
the sake of simplicity, we use the notations || • || and || • ||mo to denote the norms || • \\L'^(dvdx) 
and II • \ \L'^[Modvdx)-: respectively. Note that the two norms are equivalent since Mq is bounded 
from up and below. Similarly, the notations ((•, •)) and , •))^^ are adapted to stand for the 

inner products (•, •)L2(dt,dx) and {■,■)L^{Modvdx)^ respectively. 

We will use frequently the following facts on the sphere, for a constant vector V E M”, 


(4.1) 

(4.2) 


Vy{vV) = {Id-v(^v)V = P^±V, 
Vy ■ {{Id — V 0 v)V) = —{n — l)v ■ V, 


where V^, and V^- are used to denote the tangential gradient (and, divergence) operator on 
the sphere. Moreover, we have some other useful formulas, 


(4.3) 

(4.4) 


I Vyfdv = -{n - 1) y vfdv, 
vVy ■ A{v) dv = — J A{v) du, 


with A{v) any smooth tangent vector field. 
Now define 


{{g}) = J J 9 'ivdx, {g) = Jgdv; 
{{g))Mo = j j gM^dvdx, {g)Mo = j gModv. 


We give a lemma about the Poincare inequality on the sphere. 


Lemma 4.1 (Appendix of ([5])). We have the following weighted Poincare inequality, for 

5EF1(S-1), 

(4.5) {\^vgf)Mo > A{{g - {g)Mo f)Mo, 

where A is the Poincare constant independent ofQo. 


Due to the compatibility condition of the Hilbert expansion jfs = jf^, Hfe = we have 

(4.6) ill)Mo = ifi) = 0, {vfi)Mo = (vfi) = 0; 

(4.7) (7|)mo = (/I) = 0, (u/Dmo = (vfl) = 0. 

Hence it is convenient to introduce the mean free spaces L|^^(S"'“^) C L^(S"'“^) composed of 
functions g satisfying {g)Mo = 0- 

As pointed out in [4], the linearized operator Cq is a self-adjoint operator under the scalar 
product {91,92)12(Mo dv) = I gmModv, since we have 

{91,^092)l 2 {Mo dv) = {'^'n 9 l,'^vg 2 )L 2 (Modv) ' 

We quote the definition of an operator Cq and an equivalent Sobolev norm on the sphere 
by spectral decomposition, as in [1], 

(4^-^) Il5ll^|^^(gn-1) = {9^^o9)L2{Modv)- 

Then we have the following lemma. 




N. JIANG, L.-J. XIONG, AND T.-F. ZHANG 


Lemma 4.2 (Lemma 2 of [4]). For g € n and s > 0, we 

have 




llffll 




!)■ 


Here, we use the notation A ^ B to denote there exists a universal constant C > 0, such that 
C-^B < Al < CB. 

Furthermore, given g G there exists a constant C such that. 


(4.9) 


’^vQ)L'^(M o dv) 


^y 


5. Estimates of fi 
Recalling the equation (I2.22h of /i, 

= 9tfo + V ■ Vxfo + f?o • (P^xDo/o) , 
we introduce a new function fi satisfying /i = /iMq, then the formula 

. (M„V.A) = -Cj, 

gives the equation (the symbol ~ has been dropped for brevity): 


(5.1) 


Cofi = -^[dtfo + V • V,/o + r?o V, • (P,xZlo/o)]. 

Mq 


From the self-adjoint property of Cq, the existence and uniqueness of fi can be easily 
established by combining the Lax-Milgram theorem and the Poincare inequality, we omit the 
details. Next we give a lemma stating the boundedness of the function /i. 

Lemma 5.1. There exists a constant C > 0 depending on ||po||L°°(o,r;H'™(T)) ^.nd ||llo||L°°(o,r;H'™(T)) 
with m > 13, such that for any t G [0, T], 


(5.2) 


II/iIIfi 3 {t) + \\dth\\H^ it)<c. 

Mq Mq 


Mo 


by the Poincare inequality on the 


Proof. We split the proof into four steps. 

Step I: Estimation for ||V^/i||^|^ (hencely for 
sphere). 

To deal with the first two terms on the right hand side of the equation (15.Ij) . it follows 
from the Poincare inequality on the sphere that, 

{{Mff^dtfo + V ■ V„/o) , /l))^ = {dtfo + V ■ V:,fQ, fl)L^(dydx) 


< 


<\\dtfo + v - Vo^foW 


2 


As for the third term, we have 


go • (P^rDofo ), < r?o • (P,xPo/o)|| 

Vo ■ (-Pi;-L-^o/o)|Imo II^'!^/iIImo 


I Mo 


Mo 


< 

/'N-' 

IIVx/lllLo + hoC\\V, • (P,xPo/o)f. 

Recall that the linearized operator Cq is nonnegative, i.e., 

((£o/l , /l»M„ = ||V./i||i,„ , 
then combining the above three estimates, we get 


(5.3) 


IVo/lll 


^ < 
Mq rsj 


%fo + V ■ V,,/of +r?o||V. • (P,xPo/o)f. 
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The fact /o = PqMq yields that 

dtfo =MQdtpo + PqMov ■ dt^o, 

^xfo — xPO T PO^O^ ■ ^ 1^0) 

together with the assumption for the (SOH) model (po)^o) € L°°{0,T; H^(T)) with m > 4, 
the above two equalities result in 

\\dtfo + v-VJof <C, 

where the constant C depends on the value of ||/5 o||l°°(o,T;/I"*(t)) and ||f^o||L°°(o,r;_fi'"(T))- 
Similar argument applied to the term • (PyxL>o/o) = — (n — l)-Do^^/o + -DqV^/o shows 

m\Nv{Px±Dofo)f <c. 

Thus we are lead to the conclusion 

(5.4) II/iIImo^I|V./iIILo^^- 

Step II: Applying the operator V^; to the equation (15.Ih . and taking L^(Modudx) scalar 
product with Vxfi, we can get the estimation for (and hence, for ||Va;/i||;;^^). 

We only collect the controls for both sides of the equation (j5.1|) as follows, 

((V,£„/i, 

= ((£„V,/i, 

> ||Vt,Va;/l||^^ — II Vxfloll^oo ||V„/l||j\^^ II^i^/iIImo 
l|v.v,,/i||^^ - C liv.flollioo ||V,/i||^^, 

{{Vx{M-\dtfo + V ■ Vxfo + % V, • (Px^Dofo)]} , 

< \\Vx{M^^[dtfo + V ■ Vxfo + r/o V, • (P,xZ?o/o)]}||mo W^MImo 

<\ l|V,V,/i||^,^ + C \\Vx{M^^[dtfo + V • V,/o + • (P,xDo/o)]}||Lo • 

The above two inequalities together implies 

(5.5) - W'^v'^xflWMo ~ ll^xf^ollioo ||V.„/i||^^ 

+ \\Vx{M^^[dth + V • Vxfo + r/oV. • {P,xDofo)]}\\lj^ , 

combining with the fact ||Vt,/i||^^ < C from step I and the assumption of (po,f^o) £ 
L°°{0,T] H'^{T)) with m > 4, this gives the result 

(5.6) ||V,/i||^^<||V,V,/if^^<C. 

Step III: The estimations for ||V„clt/i||j^^^ (and ||9t/i||A,f(j) is similar to that of step II, and 
we finally get 

(5.7) ||5Ji||^^ < WVMWl,^ 

< WmoWl^ W^vfiWM, + \\dt{M^^[dtfo + V . Vxfo + . (P,xZ)o/o)]}||Lo 

<C1, 

where the constant C depends on ||po||l°°(o,T;LI"*(t)) ^^^d ||fIo||i,°°(o,T;Li"*(T)) with m > 4. 
Step IV: Performing a similar scheme, we can get the boundedness of ||V.„Va;Clt/i||^^ (and 
||Va;i9t/i||^jj) with m > 7, which yields that, together with the previous estimates, 

(5.8) II.CiIIhi^ + 11 ^ 4 / 11117 ^^ ^ c*. 

Mq Mq 

Proving the lemma for the estimates in spaces (with m > 10) and (with m > 13) 
proceeds similarly as above, here we omit the details. □ 
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Remark 5.1. Actually, we can get the similar results for higher order derivatives of fi, 

(5-9) sup iWfiWn^ (t) + WdtfiWn^ W) < C, 

te[o,r] "0 *^0 

where N >2 and C depends on H^{T) (m > AN) norm of {po, Oq). 


6. Energy Estimates for Remainder Equation 

The proof of Theorem l3.2l relies crucially on the following a priori estimate for the remainder 
equation (I2.25|) . 

(6.1) dtief) + V • V,(£/) + r?oV, • {P,vDo{ef)) + Cof = h + sfho, 

We first introduce some instant energy functionals, 

(6.2) T = To + Wi + T 2 , 0 = ^0+ ^1 + ^2) T-L = T-Lo + T~li+T~L2, 

where the terms Fi,Qi,'Hi (z = 0,1, 2) are defined as follows, 

(6.3) To = e||/||^^, eo = ||V,/||^^, To = ||/ii||Lo> 

(6.4) Ti = II + e II = £ II V,V,/||1,^ + || , 

= £ ||Va;/ji||^^ + II II , 

(6.5) F 2 = £^ Plff^, + IIV.V./IIL 3 + e \Kf\\l ^, 

^2 = II V.V2/||^^ + e ||vX/||Lo + II^'/IImo > 

T2 = II^x^iIImo + ^II^'»^*^iIIm-o + II^^^iIImo ■ 

We shall give a uniform (in e) estimate for the instant energy £, which plays the most 
important role in the process of Hilbert expansion. 


Lemma 6.1. As for the above instant energy functional, there exists an eq > 0 such that, 
for any e G (0,eo), 

(6.6) ^^£ + Q<C{£ + n), 

where the constant C depends on the value 0/||po||l°°(o,T;/I"'(t)) II^o||l°°(o,T;/I’"(t)) (with 

m > 13j. 


Proof of Lemma fTTl we shall give three type energy estimates in what follows: 

^Mo estimates, corresponding to estimates for energy functionals Fo,F 2 ,F 2 , respectively, 
and then close them to finish the proof. We also give two remarks, one of which concerns the 
estimates (for large integer N) and the other discusses the role of LL. 


6.1. Claim 1: {L'j^^ estimates). Taking L^(Modudx) inner product with /, we get from the 
equation (|6.1I1 that 


with 


|/p5fModu dx 


The second term on the left-hand side can be controllec 

1 


((t-V,/,/))mo = 


|/|^5tHo • vMq dv dx 
u 

V ■ Va;(/^Mo) du dx — 


< ||clfHo||L2 


2 

Mo¬ 


by 


v ■ VxMqI/P du dx 


1 

=»-2 


v ®v \ VxHoI/I^Mq dx dx, 


e{{vy.f,f))M,<e\\NML^\\f\\m- 


hencely 
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As for the third term, it follows that 

-Voii^v • {P,^Dof)J))Mo =%((^o/, V,/))mo =Vo [[ Do^V.fModvdx 


= -r/o(n- 1) 


9 1 

Dovf^Modvdx - -% 


Dof'^VyModv dx, 


combining with the facts that = P^±Q,qMo and |P^xflo| < 1, this gives 

|%e((V. • (P,xDo/),/))moI < ^II^oIIl- 11 / 111 , 0 . 

On the other hand, the nonnegativity of the linearized operator Cq implies that 

The Poincare inequality and the Holder inequality gives the following estimate, 

((/T,/))mo < WhiUMlMo < ||/il||Mol|V./||Mo < + ^||V,/||1,„. 

Notice that the last term on the right hand of equation ()6.1I) can be bounded as follows, 

£{{hof,f))Mo < II^oIImo ell/lllfo ^ +%Po||L-)e||/|lMo> 

then all these above estimates together gives that 
(6.7) 


1 d 


2 

Mo 


1 


+ aII^ij/IImo ~ (l|5tf^o||Lg° + ||V2;Ho||Lg° + l/o||.C)o||L^)e||/||Mo + II^iIImo- 


2 dt 
Denote 

Co = ||<9t^o||Lg° + ||Va;Do||Lg° + ) 

then there exists a universal constant C such that 


( 6 . 8 ) 


^^.^o + ^eo<C(CoT-o + 7^o). 


6.2. Claim 2: estimates). 

Step I: Taking L‘^{Mq dxdx) scalar product with To/ to control the quantity e ||V.i;/||l^^, we 
have the following estimates for the first two terms on the left-hand side, 

e{{dtf,£of))M,=e{{dtV.f,V.f))M^ 

dtMo\Vyf\'^dv dx 


_1 d 2 1 


= 7;-^^\\'^vf\\M„-l;£: II V ■ dtO .0 Mo\Vyf\'^ dvdx, 


2 dC " ' ' 2 

e ((^ • V./ , To/))mo =" ’ ^-/))Ma 

=e ((u • VJ , + e {{VJ , 


1 

~2^ 


v-\7^i\\7^f\‘^Mo)dvdx- -e // x • V^Mo| V^/p dx dx 


+ e((V./, V./)>^^ 

=0- //Mo|V^/pdxdx-be||Va:/||A^p l|V^/|| 


<e II||V„/|lt„ + Ce^- \N„f\\l,, + j ||, 

where we have used the Holder inequality and the Poincare inequality on the sphere in the 
last line. 

As for the third term, we write 

((V^ • (PyxDo/), Pof))Mo = ^^0 “ ^)Dovf + DoVyf , Tq/))^^. 


12 


N. JIANG, L.-J. XIONG, AND T.-F. ZHANG 


It is an easy matter to calculate, 

(^^/)-^o/) 1,2 (Modi;) = i'^'>}{vf),'Vvf)L2(^Modv) = if + ^’'/)l2(Mo di;) 

=\J{Vvf)Modv + I v\VJ\^Modv 

=—^ j vf^Modv - ^ j VvMof 'i^ + y v\S/vf\^MQdv. 

By Lemma 14.21 we get the fact 

(6-9) i^of ,'^vf)L^(^Modv) — ^ ll^'»/llL2(Modi;) • 

Combining with the Poincare inequality on the sphere, the two above inequalities yield 
that 

\em{{^v{P.^Dof),Cof))^J 

<em PolLgo ||V./||^^ + \\Doh^ WfWl^ 

+ er]o^ff-^ ||L>o|Iloo ||/||mo + ^Voin - 1) ||L>o|Iloo ||V^/||^^ 

<erjo II-DoIIloo ||V^/||^^ . 

On the other hand, we can control the terms on the right-hand side of equation ()6.ip as 
follows, 

((A/. £»/»„.=II v;/iiL. 

{{h , C„f))„^ = < ||V„A,||„„ ||V„/||„, < C||V„l.i|lI,„ + i IIV^/IlL . 

£ ((!.„/, C„f))„^ =e ((V./fc„, V„/»„^ + £ ((/V„fc„, 

<e II^oIIl-„ IIVi,/||mo + ^ II/IImo W^vIWmo ’ 

where we have used the Poincare inequality again. 

Accordingly, we have 

(6.10) lA(e||v„/|£,„) + l||v;/||L 

< II Vt, hi 11^^ -|- e II V3;/||j^,^^ + (||9fOo|| II V3;fIo|| + VO 11 -^0 11 ^ II^h/IImO 

+ e II^oIIlj^^ II^h/IImo + ^ l|v,;ho||ioo^ II/IImo • 

Step II: Apply the operator V^, to the equation (16.ip . and take LP‘{Mq dn dx) scalar product 
with sVxf, then we get 

e^dt^xf, W^.fWh, jj dtMo\Vxf?dvdx, 

{{y.{v • Vxf ), yxf))M, 11^- ^x{\Vxf\^Mo) dvdx 

- n • V,rMo|Va;/pdndx 

= -^e^yy v®v:Vx^oMo\Vxf\‘^dvdx. 

The third term can be decomposed into two terms, 

- e% {{VxVx • (P^xBof) , Vxf))^^ 

= - ((V. • iP^^DoVxf ), ((V. • {P^^VxDof ), V,/)>^^ 

=e^r/o(/i + h)- 
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Some straightforward calculations enable us to get 

h =((I?oV./,V,(V,/Mo))) 

={{DoVJ, ^vVxfMo + V,/V,Mo)) 


1 
2 
1 

'2 
n — 1 


Do{Vy\Vxf\‘^) Modvdx + 
Do{'^v\'^xf\‘^) Modvdx + 


Dov\VxffModvdx + - 


-Do|Va;/pV„Modn dx 
Do\Vxf\^MoP^±no dn dx 
-DolVa^/pMoP^x^o dn dx 


— 2 II-^oIIlj° II^x/IImo ’ 

and by the Poincare inequality, 

h = {{{n - l)V^Dovf - VxDoS/xf , V,/)>^^ 

<(n — 1) IIVajDollioo II/IImo II^x/IImo 11 ^ 2 ^-^ 011^00 ll^'i^/llMo II^x/IImo 
< n 11 Vj;-Do 11 2,00 11^2^/IImo ll^a^/llMo ’ 
hence the third term can be bounded by 

kNo ((V,V, • {P^xDof) , V,/))^J = e%oki + h\ 

I/O 11^0Ikoo IIVa;/!!^^ + e ho ll^x-Dolk^ ||V^,/||j^^^ HV^/IIjv^q • 

Next we turn to consider the terms on the right-hand side of equation (I6.ip . Recalling the 
nonnegativity property of the linearized operator Cq = — (A^ -|- Dq • V^), we have 

((V,£o/, eV,/))^^ =e{{CoVxf, - e((V.DoV./, V./)>^^ 

>e IIV^Vj;/!!— £ ||V2;Do|k^ II^*^/IImo II^x/IImq 


“4^ 11^ xj 


l^^-CellV.D "2 


iQII2,00 IIII 
2 


2 

Mo ’ 


where we have used the Poincare inequality on the sphere ||Va;/|||^^ < l|V.V,,/||^^. 
Note the Poincare inequality also yield that 

((Vxhi, ^’^xf ))II Va;/ii 11^^ II Va ;/||< Ce || Va;hi 11^^ -|- —e || , 

e' «V.(ho/), =£2 ((V./ho , + £2 ((/V.ho , 


<£^ 


-L?', II^2;/|Imo + ^ ll^x^olli; 


Mo 


|Vx/|| 


Mo ’ 


which ensure that, together with the previous estimates 
(d-ll) o~V^ \\''xJ\\Mo 


o-iT(e'l|Vx/|lL) + xe||V,V,/f 


Mo 

2 


2 dt'" " ■ "“0^ ' 2' 

<£ ||V3;Do|l2oo ||V^/||^^ -TO II x'nu Mo 

+ (||(9iDo|k^ -t" ||Va;Do|koo -t- i/o ll-^olk^) ^ II^x/IIa^o 

+ £ ho II Vx-DqII 2^ ll^r/llMo II^2 ;/|Ia2o 
+ (£' ll^oll2g^„ INxfWl, + £' l|v.ho||2^„ II/IImo IIVx/IUo)- 
Step III: Combining the two inequalities (Ib.lUp and (|6.1ip gives that 

k " ii^-ziIm.) +\{^ II v„v./iiC + II vhllLo) 

II^x/IImo £ II^2j/IImo) “I" II^^^iIImo) 


< 

r\j 
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+ (liatfiolLoo + llv.fiolL^ + m \\Doh^) + s II 

+ r?o||V.I?olLoo e'||V./||^J|V./||^^ 

+ {ll/^ollio.^ (e' llv./llLo + e l|V„/||Lo) + l|V./iolL» II/IImo I|V./|Imo 
+ e||/IUol|V./|lMo}- 

Denote 

Cl — sup I^Cq , II V 2 ;/l 0 II , II II ’ ^0 II VxDq II I" , 

with Cq = IIStDoll^^oo + ||Va;Do||^oo +i?o ll-Do||^oo defined as before, then we can conclude that, 
up to a constant C > 0, 

(6.13) \^^j:^ + \g^<cCi{l + e"2){Fi+Fo) + Cni, 

where it should be pointed out that 

llholli,^^ <Cq, 

I|V./||mo I|V./||mo l|V./||Mo)(e||V./||Mo) < 

6.3. Claim 3: estimates). 

Step I: Apply the operator to the equation ()6.ip . and take L^(Modndx) scalar product 
with the quantity e^V^/, then we can get immediately, 

£3 {{d,vlf, \Kf\\l^ ~¥' II |V-/l'dudx, 

{{Vliv • V,/) , JJv V,(|V2/|2)Mo dudx 

=0-^e^yy’ n-V,Mo \Vlf\^dvdx. 

Similar as before, we split the third term into two terms, 

= - e\o ((V. • [P^xiVlDof + 2V,,DoV,/)], V^/))^^ 
-e^rjo{{V,-{P,^DoVlf),Vlf))^^ 

=e^7Jo{Ih + Il2). 

It follows by calculation that, 

Ih ={{{n- IWlDovf - VlDoVJ , 

+ 2 {{{n - l)Va;DonVa;/ - V^DoVyV,;f , ^If)) 

< ((^ -1) II/IImo + I|V./|Imo) I|v^/||mo 

+ 2 llV.Dolli^^ ((n - 1) II V,/|U^ + ||V,V,/||^J || 

< II W^v/Wmo KfWuo + llVxDolLgo IIV.V./IU^ ||V^/||^„, 
Ih={{DoVlf,VyVlf))^^ = \jj DoV,(|v 2/|2) Modndx 

= ^ II Dov |V2/|2 Modvdx-^ll DoVyMo \Vlf\^dvdx 
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which give the estimate 


<e^rio \\Doh^ Kf\\l,+e^Vo I|V.V./||^^ ||v2/||^^ 

+ eSo\\vlDo\\^^\\V,fUjvlf\\^^. 

As for the right-hand side of the equation (|6.1I) . we collect their estimates, by the non¬ 
negativity property of the linearized operator Cq and the Poincare inequality on the sphere 
II^x/IImo ~ a® follows, 


{{VlCof, e^Vlf)) 

=e^ {{CoVlf , ((V^f^oV,/ + , V^/))^^ 

>^2 ||V,V2/||^^ - [ ||v2fio|L^ ||V,/||^^ PIIWm, 

+ 2 llV.Oolli^^ ||V,V,,/||^J|V2/||^^^ ■ 

>le^ ||V.V2/||^^ - Ce^ (llV^rioll^^^ W^v/Wh, + l|V.V,,/f^^ 

{{Vlh, , Vlf))^^ 


Mo 


|V.V2/|| 


Mo 


1 


72 fir 
IImo 


I Mo 

IlfcolliK. II II II V,/||„, II V=/|| 


Mo 


+ e V^j/io 


Mo 


V 


I Mo 


Then it follows from the above estimates that, 

(6.14) ~e=||Vl/||l,, + ie^||V„vr||L, 

<e"(||v 1S1„||^„ II V„/||i,„ + II V.S1„|||„ II vwnii,.) + ire'll lll,„ 

+ (llafiolli- + l|V Alli» + m l|Oolli») II Vl/lll,, 

+ e\, ||V,C„||i„ ||V„V./II„. II V2/||„^ + A,„ II V^Doll,,^ II V„/||„„ II vni„. 

+ ( llfcollig. =" II vi/HL. + II V./.olli« I|V,/||„„ II v^/||„^ 

+ ||v?l.„||,„ E^'II/IUJIVJ/II^J. 

Step II: To control the quantity ||V.i;V 3 ;/||^jj, we apply the operator Va, to the equation 
m, and take L^(Modudx) scalar product with eCoVxf, then it follows 

emdtVxf,CoVxf))^^ 

=e^ {{dtV.Vxf , V.V,/))^^ 

=\-^/ l|V.V./||L„ ~¥' II 

s\{vVlf,CoVxf))^,^ 

=e2((v,(u-V2/), V.Va./))^^ 
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+ ii^lf. v.v./))„^ 

<0 + if" ||V,f!„||i„ II V„V,/||5,, + C,’> II + ie II v;v./||"„^ , 

where we have used the Poincare inequality on the sphere ||V^V 3 ;/||jy^^ < 

The third term can be estimated by a similar decomposition as before, 

- ((V.V, • {P^^Dof), 

= - «V. • {P,±DoVJ + P,±VM), 

=e^r/o (((n - l)Dov'V^f - DoVyV^f , CqV xf)) 

+ {{{n - l)VxDovf - VxDoVyf , CqV xf)) 

=e\o{Ih + Ih). 

Combining Lemma 14.21 and the Poincare inequality on the sphere ensures that 
Ih =(n - 1) {{DoVy{v • Vxf) , V,V,/))^^ - {{DoVyVxf , 

<(n - 1) polLoo (||V,,/||^^ + ||V,V,,/||^J l|V.V,/|U„ 

+ IJ 

<{n - 1) polLoo ||V,V,/||^,^ + ll^ollioo ||V,V,/||^^ 

<n||Z)o||ioo||V,V,,/||^^, 


Ih ={n - 1) {{VxDoVyivf), VyVxf))^^ - {{VxDoVlf , 
<(n - 1) II V.Z^olLoo (II/IImo + I|V./||mo) I|V.V,/||mo 
+ \\VxDoh^ I|v^/||mo liv.v./iu^ 

<n||V.Z)olLo. ||V2/||^J|V.V./||^,^, 


therefore we get 

e^Vo ((V.V, • (Py^Dof), Co^.f))^^ 

<e% polLoo ||V,V,/||^^ + s% WVxDoh^ KAImo II^-V-ZUmo ' 


Similar discussions as before enable us to deal with the terms on the right-hand side, as 
follows. 


{{VxCof , =e {{CoVxf , - e ((V,L!oV,/ , £oV./» 


Mo 


'Mo 


l|Vi;v,/||L - e 11 V=/|| ||V„V./II„„ 


Mo 

2^I|2 


>6^l|V^V./||t,„-Ce||V.Oo||i.||v2/||^^, 

e((Va;/ii, CoVxf))j^^ =e{{VyVxhi , VyVxf))^^ < Ce ||V^Va;/ii||^Q -F 


'Mo —^ xj //^jj 

e" {{^x{hof) , ^o^J))mo {{^vVxihof), V,V,/» 


Mo 


Mo 


<e^ 


E 


|ai|-l-|a 2 |=l 
I/ 3 i| + I/ 32|=1 








Mo 


I|v.v,/|Imo • 
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All the above estimates together give that 


(6.15) 


IIVA|li„ ||v=/||i,^ + £ II 

+ (ll^t^ollLgo + ||V2;rio|lLoo + r/0 II-DqII^^oo) IIVuVoi/IImo 
+ e^o llV.Z^olLoc ||V2/||^^ ||V,V,/||^^ 


+ E 

|ai| + |a 2 | = l 

m+m=i 




v“=vfv 


I|v.v./|Imo • 

Mq 


Step III: By taking L^(Modf dx) scalar product with CqJ, we will get from the equation 
(jh.ip the control for e ||V^/||^^ in a similar process as before. Here we omit the details and 
only list the estimates, 


e{{dtf, Clf))^^ = {{dtVlf,Vlf))Mo 

= ^^e||V^/|lLo “ // 9tMo\Vlf\^dvdx, 

e {{v • V./, =e {{v • V.V^/, V^/))^^ + 2e ((V„V./, V^/))^^ 

=^e /[ ■'^x{\'^lf\^Mo)dvdx - jj • Va,Mo|V^/p du dx 

+ 2 e((V.V./,vS/»^^ 

<e II V.!iolli« II + Ce= || V.V./llL. + J II Vj/HL„, 


|e%((V,; • (P^±A)o/),iCo/))moI < {{{n - 


ll^oll 


1 ,°° 




I Mo 


l)Dovf-DoV,f,Clf))^^ 

(Using Lemma 14.21) 


((£„/, 4/))„. = II vS/llL. 

(('=1. £5 /»m. < iiv;'=iiL,. iiv;/iIm. s c iivjhinL+ \ iivj/|il , 
e ((5„/, £§/»„^ =£ ((v;(A„/), v;/»„^ 


<e E IIK''=oIIl» l|vr/l|„J|vS/||„^ 

|oi| + |«2|=2 


<E E IIK'AoIIl- 

0<|oi|<2 



2 

Mo ■ 


Combining these above estimates implies that 

(6.16) ILe ||vi/||^,^ + 1 ||V’/||L 

S||vSfti|lL + E=l|v.v./llL+ E llv?ftolL» H|v4IIL 

0<|ai|<2 

+ + llV.Llolligo + VO \\Doh^) e ||VS/||^^. 

Step IV: Summing up the inequalities (I6.14I1 . (16.151) and (|6.16p . we get 

(6.17) 

(^“ IIE/llL,+II v.v./iiL+^ II v;/|lL) 
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+ 






<E\\\vlno\\l^ W^vfWl,, + l|v,v,/||^^) + e \\vM\l^ IKfWlo 

+ + E WV^VMmo + + E^ \\Vlf\\l^ + E^ ||V,V./|| 

+ + llV.OolLg. + % \\Doh^)ie^ \Kf\\L + l|V.V,/||l,^ + e \\vlf\f ) 


2 

Mo 


+ r]o 


V.OolLj, l|V„V,/||j,„ ||v2/|| + £=> IIv^OoH,.^ iiv./iu, ||v2/|| 


Mo 


-\-E ||Va;-D( 


x-^0\\ 


Nt 


\mo ll^^’^a;/|lMo 


+ e'^ 


E 


\Pi\+m=2 


vfv 


Mo 


V 


I Mo 


+ 6 ^ E 


|ai| + |o2|=l 
I/3i| + I/32|=1 


T oo 

^x,v 




Mo 




+ e l|VrA„|li„ ||V^/|| 

0 <|«i |<2 

Denote 


2 

Mo 


C '2 = sup|c'o, ||V^D 


2 IIv 72 n ||2 

0 11 /■ oo ? ^0 /■ oo ? 


|a| + |/3|<2 


then the above equation can be rewritten as 

(6.18) 22^^ <C'C'2(1 + e + r/oe^ + r?o£)-^2 + C'H2 

+ CC2(1 + E^)F^F"i + CC2(1 + E^+ e)F^f"^ 
<C*(72(1 + 6 + rjoE^ + t]qe)F2 + CC*2(1 + e2)J^]^ 

+ C'C72(1 + e)-^o + Cil^2 
<2CC2{F2 +Fi+Fo) + Cn2. 


6.4. Closing the Estimates. Noticing the definition oi £ = Fq + Fi + F 2 , we collect the 
inequalities (j 6 . 8 p . (I6.13p . and (16.181) to get that 

(6.19) -^£ + g <C{£ + n), 

at 

holds for any e G (0,eo) and t G [ 0 ,T]. Obviously, the value of C depends upon the value 
of supjg[o^'r]{C'o, Cl, ( 72 }, which is depending eventually upon the value of \\po\\l°°{o,t-,h^{t)) 
and ||Do||ioo(Q (with m > 13). This completes the whole proof of Lemma [ 6 Tl □ 


Remark 6.1. Supposed that the initial data (pq^jDq*) are sufficiently smooth, then we can 
consider high order energy estimates for the remainder equation dEI]). Here we only give 
the result about the energy estimates because the proof is similar as before. Specifically 
speaking, 


( 6 . 20 ) 


iA 

2 dt 


^ el/3| + l 

|a|+|/3|=Af 

0<\I3\<N 



2 

Mo 



E 

\a\+\0\=N 

0<DI<1V 


V^'vf/ 


2 

Mo 


< 


EEC, 

|«| + |/3|=IV/3i+/32=/3 
i<l,a|<Ar |,ai|>i 




V^Dr 


L°° 




M 


Mo 










































+ E - 

|a| + |/3|=iV 
0<|/3|<Af 


l/3| 
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.1/31+2 


l«l + l/3|=/V 
0<|/3|<A-1 


vr'vf+v 


Mo 


Mo ^ 

+ (||aiJ^olLo= + ||V,l^o|lLo=+%PolLo.) y V“Vf/ 

l«l + l/3|=/V 
0<|/3|<A 

+ 3?o y y V“Vf/ 

|„| + |/3|=7V/3i+/3.=/3 
l<|/3|<Af |/3i|>l 

+ E E C‘c|' 

U|-|-|/3|=7V «i+a2=a 
|„i|>l /3 i+/32=/3 

+ y Cl^ 


2 

Mo 


.1/31+1 


Mo 



2 



V“vf/ 


T 00 
^x,v 


Mo 



.1/31+1 


Mo 


m=N 

/3i+/32=/3 


vfv Vf/ 

Mo 


-A +1 


Mo 


Denote 

j^N= y 

|«| + |/3|=A 
0<|/3|<Af 

= y f 

l«| + |/3|=A 
0<|/3|<A 

Cat = sup I y 


.l/3|+l 


V“Vf/ = y vr^vf/ 

Mo 




Mo 


|a| + |/3|=A 
0<|/9|<A 


2 

Mo 


C: 


'/3i 


/9i+/32=/3 

l/3i|>l 




. -loE'^s 


'/3i 

r CO ' ’ “ Z—^ /3 

" /3i+/32=/3 

l/3i|>l 


V^E»r 


. Ecj'c^' 


ai+Q;2=C»: 

/3i+/32=/3 




}. 


i/ien follows from the above inequality that 


(6.21) i <Cn{NTn + {N- l)eTN-i + • • • + 2e^-^T2 + 

+ ZZat + J-]\f + CqJ-n 

1 1 1 iV—1 AT 

+ i?oC*A-E^(A^e 2 + (.^ ~ -|-•••-|-e 2 J-i + e 2 J^g) 

+ -Eat-i + • • • + -El + -Eg) + CT-Ljq. 

Remark 6.2. In fact, recall the definition 

^ • (^.^^ 0 / 1 )] 

Mg 

= - ^ [5t/i + V ■ Vxfi - m {{n - l)Dovfi - DoVyfi)] , 

Mg 

then by simple calculations, we can infer from Lemma 15.11 that 

^o(t) = \\hi\?Mo ^ C{\\dth\\Mo > I|Vx/i|Imo ’ I|V./i|Imo) 

<C'(II^'o||l°°(0,T;A’"(T))) II^o||l°°(0,T;/R"(T)))) 

'Hl{t) =e\\S7xhl\\\i^ + IIVi^^lllLo - ^(II/^o||l°°(0,T;A"*(T))) II^^o||l°°(0,T;A’"(T)))) 

'H2{t) =e^ ||v25i||^^ +e||V,V,5i||^^ + 

<C'(II/1o||l°°(0,T;A’"(T))) II^o||l°°(0,T;A’"(T)))) 
for any t G [0, T] and m > 13. 
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As a consequence, by choosing a new constant C = supjg[o,T]{C'o) <^ 2 ) "H}, we can fet a 
refined formulation of dSSl), 

( 6 . 22 ) -^^£ + g<CC{£ + l), 

which can lead to the same result. Here we keep the original formulation with the term H in 
order to make the contributions of fi more clear. 


7. Completion of the Proof of the Main Theorem 

7.1. Local Existence of the Remainder Equation. Based on the a priori estimates 
fLemma 16.ip . we can get the local existence for the remainder equation (16.ip by a standard 
iteration scheme. 


Lemma 7.1. Given the initial datum /(O, x, v) G H‘^{T x §£ then there exist T* > 0 such 
that the remainder equation mi admits a unique solution f(t,x,v) on CPPO, T*); R^(T x 
and moreover, there exists a constant E > 0, 


sup £{f{t))<E, 
te[o,T*] 


provided that the initial datum satisfies £{f{0)) < y. 


Observing that the equation (16.Ih is linear with respect to /, proving the lemma can be 
established by a standard iteration scheme and a straightforward compactness justification, 
hence we will only sketch the proof. Here we consider the following iteration scheme: 


<9t(er+i) + V • V,,(er+i) + % V, • Do{e + £ 0 /"+^ = hi + epho, 
r+^{0,x,v) = f{0,x,v), 


starting with f^{0,x,v) = f{0,x,v). 

We remark that under the hypothesis of Theorem 13.21 there exists eo > 0 such that, 

£{P+\0)) = + E^+\0) + < I 

holds for the above E > 0 and all e G (0,eo)- 

To complete the proof of Lemma l7.ll it suffices to get a uniform (in n) estimate for £"‘^^{t). 


Lemma 7.2. There exists T* > 0, such that */sup^gjo^T*] ^(/”'(^)) ^ then 

sup £{m\t))<E. 
tG[o,r*] 

Proof. Following exactly the same lines as the proof of the a priori estimates, we can get 
finally 

< Co(T^+' + PS) + ^ 0 , 

<Ci{Pf+^+Pf + PS) + Hi, 

' < C2{PS+^ +PS + PS + PS) + ^ 2 , 

hence a constant C > 0 depending on sup^gjQ Ci, 6 * 2 } exists such that. 



(7.1) 
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Integrating from 0 to t with respect to the time variable, we have 

rt 


£{r^\t))+ [ g{r+\s))ds 
Jo 


<£{r+\ 0 )) + Ct{snp £{r+\s))+ sup £{ns)) + sup n{s)). 

sS[0,4] sG[0,t] sS[0,4] 

Take sufficiently small T* such that, 

CT* < -, Cr* sup His) < -E, 

O sG[ 0 ,t] o 

which together with the assumptions supjg[o,T*] ^(/"'(^)) ^ ^ enables 

us to get the desired result 


(7.2) 


sup £{r+\t))<E. 
te[o,r*] 


□ 

Then we are left to prove the convergence. Set tc"' = — /”, for which the iteration 

scheme dP^D leads to 

( dt{ew'^) + V • Vxiew'^) + % • {Py±Do{sw'^)) + Cqw^ = ew^~^ho, 

\ w'^{0, X, v) = 0. 


(7.3) 


The same computations as used for ( |PnD give that 

£{w'^{t)) + f G{w^{s)) ds < CT^,{ sup £{w"'{s)) + sup £(s))). 

Jo sG[0,T*] 

Take T* small enough so that CT* < i, then we obtain 

1 


sG [OjT** 


sup £{r+\t) - nt)) < - sup £{r{t)-r-\t)). 
te[o,n] ^ te[o,r*] 


Consequently, it can be proved easily the sequence {/”} is convergent in the space C'(0, T*; 
77^(T X S”"”^)) and so is its limit /. This completes the proof of Lemma l7.11 

7.2. Completion for the Main Theorem. We are now in a position to complete the proof 
of our main Theorem El Lemma rm shows that there exists a maximal time > 0 such 
that the solution /| to the remainder equation (|6.1I) satisfies the a priori estimate stated in 
lemma 16.11 

(7.4) ^e^Q<c{£p'K), 

for any e G (0,eo) and t € [0,T£]. 

Observing the hypothesis of Theorem 13.21 11 /111^2^(0) < C implies that T(0) < E. Define 

E 2 = e^'^(P(0) + CT sup P(t)), 
te[o,T] 

To = sup{t G [o,r£], £{t) <E 2 }. 


We claim that > T. Indeed, applying the Gronwall inequality, we infer from the above 
a priori estimate, for t G [0,T£], 

(7.5) £{t) < e^*(P(0) + Ct sup P(s)). 

If Tg < T, then it follows that 


£{t) < E2, 
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which in turn yields that = Tq, and furthermore, the a priori estimate are adapted to the 
time t = Tg so that the solution may be continued beyond T^. This stands in contradiction 
to the maximal property of Tg. Thus Tg > T, and £{t) < E 2 holds for t G [0,T]. 

Recalling the definition of £ 


£ = J-Q + J-i + J-2 

=e 


' + (s' llV./llLo + e ||V2/||' + £2 ||V,V./||^^ + £ ||V2/|| 


Mo 


and the equivalence between the two norms 

o, 


I Mo 


V - XJ iiMo ' ^ W' VJ WMq) 

and II • II, we actually have proved the result 


(7-6) £2 11/111^2^^ + £||/|||_h-i_^ + £2 II/, 

with C is independent of £ G (0,£o) and t G [0,T]. 


2\\Hl, ^ 


< c, 
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